NEWTON POLYTOPES AND WITNESS SETS 
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Abstract. We present two algorithms that compute the Newton polytope of a polyno- 
mial defining a hypersurface % in C" using numerical computation. The first algorithm 
assumes that we may only compute values of / — this may occur if / is given as a straight- 
line program, as a determinant, or as an oracle. The second algorithm assumes that "H 
is represented numerically via a witness set. That is, it computes the Newton polytope 
of H using only the ability to compute numerical representatives of its intersections with 
lines. Such witness set representations are readily obtained when H is the image of a 
map or is a discriminant. We use the second algorithm to compute a face of the Newton 
polytope of the Liiroth invariant, as well as its restriction to that face. 



Introduction 

While a hypersurface % in C™ is always defined by the vanishing of a single polynomial 
/, we may not always have access to the monomial representation of /. This occurs, for 
example, when T-L is the image of a map or if / is represented as a straight-line program, 
and it is a well-understood and challenging problem to determine the polynomial / when % 
is represented in this way. Elimination theory gives a symbolic method based on Grobner 
bases that can determine / from a representation of Ji as the image of a map or as a 
discriminant (?]. Such computations require that the map be represented symbolically, 
and they may be infeasible for moderately-sized input. 

The set of monomials in /, or more simply the convex hull of their exponent vectors (the 
Newton polytope of /), is an important combinatorial invariant of the hypersurface. The 
Newton polytope encodes asymptotic information about H and determining it from % is a 
step towards determining the polynomial /. For example, numerical linear algebra [5} fTTj 
may be used to find / given its Newton polytope. Similarly, the Newton polytope of an 
image of a map may be computed from Newton polytopes of the polynomials defining the 
map [TOJ [121 E2J 123 [26] , and computed using tropical geometric algorithms [27] . 

We propose numerical methods to compute the Newton polytope of / in two cases 
when / is not known explicitly. We first show how to compute the Newton polytope when 
we are able to evaluate /. This occurs, for example, if / is represented as a straight- 
line program or as a determinant (neither of which we want to expand as a sum of 
monomials), or perhaps as a compiled program. For the other case, we suppose that / 
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defines a hypersurface % that is represented numerically as a witness set. Our basic idea 
is similar to ideas from tropical geometry. The tropical variety of a hypersurface "H in 
(C x ) n is the normal fan to the Newton polytope of a defining polynomial /, augmented 
with the edge lengths. The underlying fan coincides with the logarithmic limit set [3j H] 
of "H, which records the asymptotic behavior of "H in (C x ) n . We use numerical nonlinear 
algebra to study the asymptotic behavior of H in (C x ) n and use this to recover the Newton 
polytope of a defining equation of "H. As both our algorithms are easily parallelizable, this 
numerical approach to Newton polytopes should allow the computation of significantly 
larger examples than are possible with purely symbolic methods. 

This paper is organized as follows. In Section (TJ we explain symbolic and geometric- 
combinatorial preliminaries, including representations of polytopes, Newton polytopes, 
and straight-line programs. In Section El we discuss the essentials of numerical nonlinear 
algebra (also called numerical algebraic geometry [21]), in particular explaining the fun- 
damental data structure of witness sets. Our main results are in the next two sections. In 
Section [3] we explain (in Theorem H] and Remark \5§ how to compute the Newton polytope 
of /, given only that we may numerically evaluate /, and in Section HJ we explain (in The- 
orems [H] and O and Remark [TP]) how to use witness sets to compute the Newton polytope 
of /. Illustrative examples are presented in these sections. In Section [5j we combine our 
approach with other techniques in numerical nonlinear algebra to explicitly compute the 
hypersurface of even Liiroth quartics. 

1. Polynomials and Polytopes 
We explain necessary background from geometric combinatorics and algebra. 

1.1. Polytopes. A polytope P is the convex hull of finitely many points A C M n , 

(1.1) P = conv(A) :={^A Q a : \ a > , ^ X a = l} . 

Dually, a polytope is the intersection of finitely many halfspaces in IR n , 

(1.2) P = {x E R n : Wi-x < k for z = l,..., JV}, 

where Wi, . . . , Wn £ K n and bi, . . . ,bjy E.WL. These are two of the most common representa- 
tions of a polytope. The first (II. ip is the convex hull representation and the second (11.21) 
is the halfspace representation. The classical algorithm of Fourier-Motzkin elimination 
converts between these two representations. 

The afhne hull of a polytope P is the smallest affine-linear space containing P. The 
boundary of P (in its affine hull) is a union of polytopes of smaller dimension than P, 
called faces of P. A facet of P is a maximal proper face, while a vertex is a minimal face 
of P (which is necessarily a point). An edge is a 1-dimensional face. 

In addition to the two representations given above, polytopes also have a tropical rep- 
resentation, which consists of the edge lengths, together with the normal fan to the edges. 
(This normal fan encodes the edge-face incidences.) Jensen and Yu [T^] gave an algorithm 
for converting a tropical representation into a convex hull representation. 
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Every linear function i 4 w • i on K" (here, w G M n ) achieves a maximum value on 
a polytope P. The subset P w of P where this maximum value is achieved is a face of 
P, called the face exposed by w. Let hp(w) be this maximum value of w ■ x on P. The 
function w H- h P {w) is called the support function of P. The support function encodes 
the halfspace representation as 

P = {iGK" : w ■ x < h P (w) for toGl"}. 

The oracle representation is a fourth natural representation of a polytope P. There are 
two versions. For the first, given w 6 1", if the face P w exposed by w is a vertex, then it 
returns that vertex, and if P w is not a vertex, it either returns a vertex on P w or detects 
that P w is not a vertex. Alternatively, it returns the value hp{w) of the support function 
at w. The classical beneath-beyond algorithm [HJ §5.2] uses an oracle representation of 
a polytope to simultaneously construct its convex-hull and halfspace representations. It 
iteratively builds a description of the polytope, including the faces and facet-supporting 
hyperplanes, adding one vertex at a time. The software package iB4e [18] implements this 
algorithm. Another algorithm converting the oracle representation to the convex hull and 
halfspace representation is "gift-wrapping" [5]. 

Our numerical algorithms return oracle representations. 

1.2. Polynomials and their Newton polytopes. Let N = {0, 1, . . . } be the nonneg- 
ative integers and write C x for the nonzero complex numbers. Of the many ways to 
represent a polynomial / G C[xi, . . . , x n ], perhaps the most familiar is in terms of mono- 
mials. For a G N n , we have the monomial 

— ai a 2 _ a n 

which has degree \a\ := ai+- ■ - + a n . A polynomial / is a linear combination of monomials 

(1.3) / = ^c a x a c a £C, 

where only finitely many coefficients c a are nonzero. The set {a G N n : c a ^ 0} is the 
support of /, which we will write as A(f), or simply A when / is understood. 

A coarser invariant of the polynomial / is its Newton polytope, N(f). This is the 
convex hull of its support 

NU) ■= oonv(A(f)). 
For w G M", the restriction f w of / to the face Af(f) w of Af(f) exposed by w is 

(1.4) /„ := Yl c - x ^ 

the sum over all terms c a x a of / where w ■ a is maximal (and thus equal to h^(f\(w).) 

A hypersurface H C C n is defined by the vanishing of a single polynomial, % = V(/). 
This polynomial / is well-defined up to multiplication by non-zero scalars if we require it 
to be of minimal degree among all polynomials vanishing on %. We define the Newton 
polytope, Af(H), of H to be the Newton polytope of any minimal degree polynomial 
/ G C[xi, ... ,x n ] defining %. 
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Polynomials are not always given as a linear combination of monomials (jl.3p . For exam- 
ple, a polynomial may be given as a determinant whose entries are themselves polynomials. 
It may be prohibitive to expand this into a sum of monomials, but it is computationally 
efficient to evaluate the determinant. For another example, a polynomial may be given 
as an oracle or as a compiled program. 

An efficient encoding of a polynomial is as a straight line program. For a polynomial 
/ : C n -> C, this is a list 

(f-n, ■ ■ ■ , f-1, fo, fl, ■ ■ ■ , fl) 

of polynomials where f — fi and we have the initial values /_j = x» for % = 1, . . . , n, and 
for every k > 0, /& is one of 

/i + /j , /i • / 3 - , or c , 

where i,j<k and c G Q[a/— I] is a Gaussian rational number. (Gaussian rational numbers 
are used for they are representable on a computer.) 

Our goal is twofold, we present an algorithm to compute the Newton polytope of a poly- 
nomial / that we can only evaluate numerically, and we present an algorithm to recover 
the Newton polytope of a polynomial / defining a hypersurface % that is represented 
numerically as a witness set (defined in §[2] below). 

In the first case, we explain how to compute the support function hj^j-in of the Newton 
polytope of /, and to compute N{f) w , when this is a vertex. This becomes an algorithm, 
at least for general w, when we have additional information about /, such as a finite 
superset B C Z n of its support and bounds on the magnitudes of its coefficients. This is 
discussed in Remark [51 

In the second case, we show how to compute Af(f) w , when this is a vertex. This is 
discussed in Remark [TUl 

2. Numerical nonlinear algebra and witness sets 

Numerical nonlinear algebra (also called numerical algebraic geometry [21]) provides 
methods based on numerical continuation for studying algebraic varieties on a computer. 
The fundamental data structure in this field is a witness set, which is a geometric repre- 
sentation based on linear sections and generic points. 

Given a polynomial system F: C m — > C n , consider an irreducible component V C 
V(F) := F" 1 ^) of its zero set of dimension k and degree d. Let C: C m — > C fc be a 
system of general affine-linear polynomials so that V(C) is a general codimension k affine 
subspace of C m . Then W := V fl V(£) will consist of d distinct points, and we call the 
triple (F, £, W) (or simply W for short) a witness set for V. The set W represents a 
general linear section of V. Numerical continuation may be used to follow the points of 
W as £ (and hence V(£)) varies continuously. This allows us to sample points from V. 

Ideally, V is a genetically reduced component of the scheme V(F) in that the Jacobian 
of F at a general point w G W C V of V has a fc-dimensional null space. Otherwise 
the scheme V(F) is not reduced along V. When V is a generically reduced component 
of V(F), the points of W are nonsingular zeroes of the polynomial system [£]. When 
V(F) is not reduced along V, the points of W are singular zeroes of this system, and it is 
numerically challenging to compute such singular points. 
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The method of deflation, building from [22], can compute W when V(F) is not reduced 
along V. In particular, the strong deflation method of [17] yields a system F' : C m — > C n 
where n' > n such that V is a generically reduced component of the scheme V(F'). 
Replacing F with F', we will assume that V is a generically reduced component of V(F). 

The notion of a witness set for the image of an irreducible variety under a linear map was 
developed in [13]. Suppose that we have a polynomial system F: C m — > C n , a generically 
reduced component of V(F) of dimension and degree d, and a linear map u : C m — > C p 
defined by cu(x) = Ax for A £ C pxm . Suppose that the algebraic set Z7 = w(V) C C p has 
dimension fc' and degree d'. A witness set for the projection U requires an afline-linear 
map £ adapted to the projection u>. Let B be a matrix ] where the rows of the matrix 
.Bi £ c fc ' xm are general vectors in the row space of A and the rows of B 2 £ C ( - fc ~ fc '- )xm are 
general vectors in C m . Define £: C m -»■ C fc by £(x) = 5x - 1 and set W := V n V(£). 
Then the quadruple (i 71 , a;, £, W) is a witness set for the projection U. By our choice of 
B, the number of points in uj(W) is the degree d 1 of U and for any fixed u £ u;(W), the 



number of points in W H u" 



[u 



is the degree of the general fiber of u restricted to V. 



Note that k — k' is the dimension of the general fiber. 

Example 1. Consider the discriminant hypersurface Ji C C 3 for univariate quadratic 
polynomials, that is, % := V(/) where f(a, b, c) = b 2 — 4ac. The triple (/, £, W) where 



£(a, b, c) 



2a -26 + 3c - 1 
3a + b — 5c — 1 



and W = Tin V(£), which consists of the two points, (a, 6, c), 

{(0.3816, -0.1071, 0.00752) , (1.2243, 2.1801, 0.97058)} 
is a witness set for %. 



This discriminant also has the form % = w(V) where u is the linear projection mapping 
(a, 6, c, x) to (a, 6, c) and V = V(F) where 



F(a, b, c, x) 



ax 2 + bx + c 
2ax + b 



This variety V has dimension 2 and degree 3, and u is defined by the matrix 

A-- 



10 
10 
10 



The quadruple (F, u, £', W) where C'(a,b,c,x) = £(a,6,c) and W = V n V(£'), which 
also consists of two points, (a, 6, c, x), 

{(0.3816, -0.1071, 0.00752, 0.1403882) , (1.2243, 2.1801, 0.97058, -0.8903882)} , 

is also a witness set for %. In particular, u(W') = W and we see that u restricted to V 
is generically one-to-one. □ 
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3. Newton polytopes via evaluation 

We address the problem of computing the Newton polytope of a polynomial / : C n — > C 
when we have a method to evaluate /. This is improved when we have some additional 
information about the polynomial /. 

For t a positive real number and w G M. n , set t w := (t Wl ,t W2 , . . . ,t Wn ). Consider the 
monomial expansion of the polynomial /, 

/ = c a x a where c a G C x . 

a£A 

For x G C n , we define 

t W .X '.= [t Wl X\ , t W2 X2 , . . . , t Wn X n ) , 

the coordinatewise product, and consider the evaluation, 

(3.1) f(t w .x) = J]c/V. 

aeA 

Let T := ftf(f) w be the face of Af(f) that is exposed by w. Then if a G J 7 , we have 
w ■ a = hj\fm(w). There is a positive real number d w such that if a G A \ J 7 , then 
w ■ a < hjsf(f)(w) — d w . Thus (13. ip becomes 

f(t w . x ) = ^ c a t wa x a + c a t wa x a 
aeAnJ 7 aeA-^J 7 

= t h ^ w) (f w {x) + ^2 cj w - a ~ h ^^ {w) x a 

aeA^J 7 

where f w is the restriction of / to the face J 7 . Observe that no exponent of t which occurs 
in the sum exceeds —d w . This gives an asymptotic expression for t 3> 0, 

(3.2) log|/(r.x)| = h m) (w)\og(t) + log|/ w (x)| + 0(t~ d -), 
from which we deduce the following limit. 

Lemma 2. If f w (x) ^ 0, then 

hm)iw) = i im . 

t/n 1 i->oo log(t) 

Thus we may approximate the support function of M{f) by evaluating / numerically. 

Remark 3. To turn Lemma |5] into an algorithm for computing h/j-/f\, we need more 
information about /, so that we may estimate the rate of convergence. For example, if 
we have a bound, in the form of a finite superset B C N n of A, then {w ■ a \ a G £>} is a 
discrete set which contains the value of htf(f)(w), and therefore the limit in Lemma [2j □ 

When w is generic in that a i— > w ■ a is injective on A, then the face Af(f) w of A/"(/) 
exposed by w is a vertex so that /^(x) ^ for any x G (C x ) n . We may dispense with the 
limit given an a priori estimate on the magnitude of the coefficients of /. 
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Theorem 4. Let f : C n — > C be a polynomial with monomial expansion, 



Suppose that 5, A > 1 and i3 C N n are snc/i i/jai 

(1) log |c Q | < 5 /or a// a E A, 

(2) log |c Q | < A + log jc^l /or all a, (3 <E A, and 

(3) Ac B with \B\ < oo. 
Let w G M. N be general in that 



d„ 



mm 



W 



a-w ■ f3\ > 



Then the face of M(f) exposed by w is a vertex which equals the unique (3 G B such that 



log mm 



log(t) 



< 



d 



w 

2 



where t > is any number with log(t) exceeding max{2A, 2(5 + e 1 ), A + log \B\ + l}/d v 
Similarly, the face exposed by —w is a vertex which equals the unique (3 G B such that 



-w(3 



iog|/(r 



log(t) 



< 



where t > is any number with log(t) exceeding max{2A, 2(5 



A + log|i3| + l}/d u 



Remark 5. Suppose that we know or may estimate the quantities B, 5, and A of 
Theorem HI Then, for general w G M n we may compute d w , and therefore evaluating 
log |/(f)|/log(i) for t dw > max{e 2A , e 2+2<5 , |i3|e A+1 } will yield w ■ (3 and hence (3. 

Even without this knowledge, we may still compute the support function h^^(w) for 
w G Q n as follows. For 0^w6 Q n the map Z n — > Q given by (3 H- w ■ (3 has image a free 
group Zd w for some d w > 0. For x G C n with f w (x) ^ and t := e T with r > 0, we have 



log|/(e™.x)| 



T 



h M(f)( 



W 



log |/«,(ar)| 



0(e 



Since hj^if\(w) G 1<d w , we may do the following. Pick a general x G C n (so that 
fw(x) 7^ 0), and compute the quantity 



(3.3) 



log|/(e™.x)| 



for r in some increasing sequence of positive numbers. We monitor ( 13. 3p for ^-convergence 
to some nd w G Zd w . Then h^f(f)(uj) = nd w . 
Every such computation gives a halfspace 

{x G M. n | w ■ x < htf(f){w)} 

containing jV(/). Since Af(f) lies in the positive orthant, we may repeat this one or more 
times to obtain a bounded polytope P containing Af(f). Having done so, set B := PflN n . 

Suppose that w G M. n is general in that the values of w • a for a G B are dis- 
tinct. This implies that w exposes a vertex /3 of Af(f)- Then a similar (but simpler 
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as f w {t w ) = cpt hM ^^ w ^) scheme as described above will result in the computation of the 
support function h^f)(w) and the vertex (3. □ 

Proof of Theorem® By the choice of w, the face of J\f(f) it exposes is a vertex, say (3 G A, 
and we have w ■ (3 = hj^^i^w). We may write 



Cpt 

Taking absolute value and logarithms, and using that log \cp\ < S and w ■ (3 = h(w), 



(3.4) 



log|/(t* 



loglc^l+log 



1 f(t w ) ~ opt 



< 5 + w ■ (3 log(t) + log 



cpt w -P 

1 , f(t w )-cpt 
cpt w -? 



w-0 



Let us estimate the last term. As A C £>, we have 



f(t w )-cpt w -P 



t w-p 



a +P 



< ^^e x t~ dw < \B\ e A_log W dw . 



aeA 

a^p 



Since log(t) > (A + log \B\ + l)/d w , we have \B\e x ~ log ^ dw < e~ l . Since log |1 + x\ < \x\ 
we have 



log 



f(t w ) - Cpt^ 



CpP 



< 



f(t w ) - cpt w -P 



t wp 



< \B\ e x ~ dwloe{t) < e -1 . 



Finally, as we have log(t) > 2(5 + e 1 )/d w , we obtain 
(3.5) 



log(t) log(t) 2 



For the other inequality, using (13 .4p and Condition (2) of the theorem, 

f{t w ) - cpt w -P 



(3.6) 
Since 



log 1/(01 > 5-A + log(t) W -/3 + log 

f(t w ) - cpr-e 



cpt w - 



w-P 



cpt w - 



w-P 



< e 



the logarithm on the right of (I3.6P exceeds —1. As 5 — A > 1 — A, we have 



iog|/(QI 

log(t) 



> 



- — j-r + w-13 > w-p- — 
log(t) 2 



since d w \og(t) > 2A. Combining this with (13. 5p proves the first statement about f(t w ). 
The statement about f(t~ w ) has the same proof, replacing w with —to. □ 
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Example 6. Reconsider the polynomial f(a, b, c) = b 2 — Aac from Ex. [T] with the vector 
w = (—1.2, 0.4, 3.7). Suppose that we take A = 5 = 2 and B = {a 2 , ab, ac, b 2 , be, c 2 } which 
are the columns of the matrix 

/2 1 1 N 
B=010210 
\0 1 1 2 / 

Then the dot products are w ■ B = (—2.4, —0.8,2.5,0.8,4.1,7.4), so that d w = 1.6. Since 
we need log(t) > 3.75, we can take t = 45, and so t w = (45~ L2 , 45 a4 , 45 3,7 ). We compute 

l0 f (Q! = 2.864 and -MM^} ^ . 8016 . 
log(t) log(t) 

Thus, the monomials ac and b 2 are the vertices J\f(f) w and J\f(f)- W , respectively. □ 

4. Newton polytopes via witness sets 

Let % C C n be an irreducible hypersurface and suppose that we have a witness set 
representation for H. As discussed in Section [2j this means that we may compute the 
intersections of Tini where I is a general line in C n . We explain how to use this information 
to compute an oracle representation of the Newton polytope of H. 

The hypersurface % C C n is defined by a single irreducible polynomial 

(4.1) / = J]c Q x a c tt eC x , 

which is determined by Ti up to multiplication by a scalar. 

Let a, b G C n be general points, and consider the parametrized line 

l a h = £( s ) ■= { sa - b | s G C} . 

Then the solutions to f(£(s)) = parameterize the intersection of H with the line 
which is a witness set for H. 

Let w G M n . For t a positive real number, consider f(t w .£(s)), which is 

(4.2) c -( Sfl i - b i) ai ( sa 2 - b2T 2 ■ ■ ■ (sa n - b n ) a - .t w - a 

Write (sa — b) a for the product of terms (sa, — bi) ai appearing in the sum. 

Let T := Af(H) w be the face of the Newton polytope of H exposed by w. If a G J 7 , 
then w ■ a — h(w), where h is the support function of Af(7-L). There is a positive number 
d w such that if a G A \ J 7 , then w ■ a < h(w) — d w . We may rewrite (14. 2p . 

f(t w l{s)) = t h ^ c a (as-b) a + c a {as - b) a t w ' a . 

Multiplying by t~ h ^ and rewriting using the definition ( j 1.4ft of f w gives 

(4.3) r h ^f(t w l(s)) = f w (£(s)) + c a {as - b) a t w a - h{w) . 

Observe that the exponent of t in each term of the sum over A \ J 7 is at most —d w . 
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As s 14 £(s) and s H- t w .£(s) are general parametrized lines in C n , the zeroes (in s) of 
f(t w .£(s)) and f w (£{s)) parameterize witness sets for / and respectively. The following 
summarizes this discussion. 

Lemma 7. In the limit as t — > oo, there are deg(f) — deg(f w ) points of the witness set 
f(t w .£(s)) = which diverge to oo (in s) and the remaining points converge to the witness 
set f w (t(s)) = 0. 



When N(H) W is a vertex fa then f w = cpx 13 (and deg(/ u 



fMs)) = c ( s(sa 1 -b 1 ) l ^(sa 2 -b 2 )^---( 



sa r 



, and 

cp{sa — b) 13 



In particular, there will be fa points of f(t w .£(s)) = converging to bi/ai as t — > oo, and 
so Lemma [7] gives a method to compute the vertices /3 of J\f{H). We give some definitions 
to make these notions more precise. 

Let a G (C x ) n and b G C n be general in that the univariate polynomial f{£ a ,b{s)) has 
d = deg(K) nondegenerate roots, and if i ^ j, then bi/ai ^ bj/aj. For any w G M. n with 
J\F(H) W = consider the bivariate function g a ,b,w(s,t) = g(s,t) := f (t w .£ a ^(s)) . Since 
g(s, 1) has d simple zeroes, there are at most finitely many positive numbers t for which 
g(s,t) does not have d simple zeroes. Therefore, there is a t > and d disjoint analytic 
curves s(t) G C for t > to which parameterize the zeroes of g{s,t) for t > to (that is, 
g(s(t),t) = for t> to). 

By Lemma [7] and our choice of a, b, for each i = 1, . . . , n, exactly fa of these curves will 
converge to bi/ai as t — )• oo, for each i = 1, . . . , n, while the remaining d — \(3\ curves will 
diverge to infinity. We give an estimate of the rates of these convergences/divergences. 

Let w G W 1 be general in that Af(7i) w is a vertex, fa Let d w be as above, and set 



C 



max{|c a | : a G A} 



Furthermore, set a min := min{l, |aj| : i — 1, . . . ,n}, a 
and the same, b min and b max , for b. Finally, for each i 



max 
1, 



Ji 



max 



1 


h 


_h 










a i 


' 2 


b ± 


_h 


, "max 







max{l, \ai 
,n, define 

and 



n}, 



We give two results about the rate of convergence/divergence of the analytic curves s(t) 
of zeroes of g(s,t), and then discuss how these may be used to compute M{%). 

Theorem 8. With the above definitions, suppose that s: (t ,oc) — > C is a continuous 
function such that g a ,b,w(s(t),t) = for t > t and that s(t) converges to bi/ai as t — >■ oo. 
Let t\ > to be a number such that if t > t\ then 



(4.4) 



s(t) ^ 



< 



li 
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Then, for all t > t\, 



(4.5) 



s(t) 



Cl; 



* < t~ d - -C-\A\- (— + - 

\ ^min \ li 



Theorem 9. With the above definitions, suppose that s: (to,oo) — > C is a continuous 
function such that g a ,b,w{s{t),t) = for t > to and that s(t) diverges to oo as t — >■ oo. Let 
t\ > t be a number such that ift > t\ then 

(4.6) \s(t)\ > > 2. 

<^min 

Then, for all t > t\, 

(4.7) \s{t)\ d ~ m > 



C-\A\ \2(a max + a min ) 

Remark 10. Theorems [8] and [9] lead to an algorithm to determine vertices of N(Ji). 
First, choose a, b G C" as above and compute 7,, 6 max , and a m i n - For a general w G M. n , 
follow points in the witness set H D (£™.^<i,&(s)) as i increases until the inequalities (I4.4p 
and (14. 6p are satisfied by the different points of the witness sets, at some t\. This will give 
likely values for the integer components of the vertex /3 exposed by w. Next, continue 
following these points until the subexponential convergence in ( I4.5P and (14. 7p is observed, 
which will confirm the value of /3. 

If we do not observe clustering of points of the witness set at s — b^/di and s = 00, 
then we discard w, as it is not sufficiently general. That is, either it exposes a positive 
dimensional face of A/"("H) or else it is very close to doing so in that d w is too small. □ 

Proof of TheoremU Fix t > t x . Since = g a ,b,w(s(t),t) = f(t w .l a , b (s(t))) and f w (x) = 
cpx 13 , (14. 3 p gives 

\(s(t)a-bf\ < £ t^-^- l -^\.\(s(t)a-br\ 

(4.8) < r*»-c- Yl \( s ( t ) a - b T\ ■ 

aSA\{/3} 

For any % and j we have 

Is^dj-bj] = \ aj \ ■ \s(t) - || < <w|s(t)-£ + |-|| < a ma x(7i + r 4 ). 
Since 2 < a max Fj and if a G A, then \a\ < d, we have 

(4.9) |(s(t)a-&) a | < (0^(74 + TO)*. 
With (I4.8p . this becomes 

(4.10) \(s(t)a-bf\ < r*™ ■C-\A\-(a max (li + r i )) d . 
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If j 7^ i, then 



s(t)a 3 -b 3 \ = \ aj \.\s{t)- b f\ = | fli |.| a (t)_jL + jL_5 



h h — b. 



> a r 



cbj I I v ' I 

> a min ■ (27^ - 7j) = a min 7j . 



Since a m ; n 7j < 1 and < d, we have 
(4.11) 



Observe that we have 



ft 



III " 6 i) ft | > K^)^ 



i(t)oi - 



\{s{t)a-bf\ 



Combining this with ( I4.10p and ( 14. lip gives 

bi' 



s(t) 



a; 



d 1 



Since 1 > a min > 7; and d > fa > 0, we have 



s(t) ■ ^ 



ft 



„ft f . -.W-ft 

a min l u mm /«/ 



i + — 



< r dw -c-\A\ 

\ u min \ Ji J J 

which completes the proof. □ 
Proof of Theorem^ Fix t > t\. Then \s(t)\ > 2. Since g a ,b,w(s(i),t) = 0, we have 
\c p t w - p {s{t)a - bf\ = | c a t w a {s{t)a - b) a < ^ \cj w ' a (s(t)a - b) a \ . 

a£yK{/3} aeA\{/3} 

Factoring out powers of \s(t)\, we obtain 

t w -P\s{t)\W\cp\\(a-bs(t)-y\ < * w ' a l 5 (<)l |a| l c all(a-6s(<)" 1 ) a |- 

a£A\{/3} 

Since s(t) — > oo as t — > oo, we must have d — \/3\ > 0. Dividing by most of the left hand 
side and by |s(t)| d and using the definition of d w , we obtain 



a&A\{/3} 



■0 ln ( f)l \ a \-d\ c «\ \(a-bs(t) 

Ul \ cp \ \(a-bs(tyiy\ 



(4.12) 



< f-dv.n. V \ n (f)\M-d \( a - bs ( t ) 
" ° 2^ |( a _ 6a ( t )-i)/»| 

a£A{(3} lv 

We estimate the terms in this last sum. As \s(t)\ > 2, for any z we have 

K - &i Sity 1 ] < \di\ + |6j s(t) _1 | < flmax + ^max, 
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and so \(a — bs(t) 1 ) Q | < (a max + b max )^. Similarly, for any i we have 

(Zmin ^ 



\ai-bis(t) : | > \ai\-\bis(t) x | > a min - b 



mm "mm 

min "max 



IK 



Thus 

|(a-&s(*) _1 )1 . , , ,| Q |, 2 (2{a msx + b* 



^ < (a m » + L x ) H (-)'" < (- 



\{a-bs{ty i )l 

Substituting this into (14. 12ft completes the proof of the theorem. □ 



Example 11. We demonstrate the convergence and divergence bounds by considering 
the polynomial f(x,y) = x 2 + 3x + 2y — 5 with the hypersurface "H := V{f) it defines. 
We have A = {1, x, y, x 2 } with \A\ = 4 and will take C = 5, a = (2 + 3 - 2 V /T T), 
b = (— 1— V— T, 2-3-\/ z r), a min = 1, a max = y/l3, b min = 1, and 6 max = Additionally, 
7i = I\ « 1.5342 for i = 1,2. 

First, consider the vector w = (1,1) for which M{%) w = (2,0) and d w — 1. We 
have g a ,b(s,t) = f(t ■ (sai — bi),t ■ (sa 2 — b 2 )), and g a ,b(s,t) = has two nonsingular 
solutions for all t > 0. Since Af(7i) w = (2,0) both solutions paths converge to &i/cti 
as t — y oo. The following table compares the actual values for the two solution paths, 
Si(t) and S2(t), with the upper bound (14.51) in Theorem [SJ In particular, this table shows 
\si(t) — bi/ai\ 2 ~ 2.2t~ 1 whereas the upper bound is 1040t~ 1 . 



t 


\si(t) - h/a^ 2 


\s 2 (t) - h/atl 2 


Upper bound (14.51) 


le2 


0.26 


0.19 


10.4 


le4 


2.2e-4 


2.2e-4 


0.104 


le6 


2.2e-6 


2.2e-6 


1.04e-3 


le8 


2.2e-8 


2.2e-8 


1.04e-5 



We now consider the vector w = (—1, —1) for which Af(7-L) w = (0, 0) and d w = 2. With 
the same a, b as above, g a ,b(s,t) = f(t~ l ■ (sa\ — bi),t~ l ■ (sa 2 — b 2 )) and g at b(s,t) = 
has two nonsingular solutions for all t > 0. Since N(H) W = (0,0), both solution paths 
diverge to oo as t — > oo. The following table compares the actual values for the two 
solution paths, s\(t) and s 2 (t), and the lower bound (14. 7p in Theorem |9j This table shows 
|si(£)| 2 « t 2 /8.71 whereas the lower bound is t 2 /4160. 



t 


l«i(*)| a 


\s2(t)\ 2 


Lower bound (14.71) 


le2 


1.17e3 


1.13e3 


2.40 


le4 


1.15e7 


1.15e7 


2.40e4 


le6 


1.15ell 


1.15ell 


2.40e8 


le8 


1.15el5 


1.15el5 


2.40el2 



□ 
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5. Even Luroth quartics 

Associating a plane quartic curve to a defining equation identifies the set of plane 
quartics with P 14 . This projective space has an interesting Luroth hypersurface whose 
general point is a Luroth quartic, which is a quartic that contains the ten vertices of some 
pentalateral (arrangement of five lines). The equation for this hypersurface is the Luroth 
invariant, which has degree 54 [2T] and is invariant under the induced action of PGL(3) 
on P(S* 4 C 3 ) ~ P 14 . A discussion of this remarkable hypersurface, with references, is given 
in [SI Remark 6.3.31]. 

We use the algorithm of Section H] to investigate the Luroth polytope, the Newton 
polytope of the Luroth invariant. While we are not yet able to compute the full Luroth 
polytope, we can compute some of its vertices, including all those on a particular three- 
dimensional face. This face is the Newton polytope of the Luroth hypersurface in the 
five-dimensional family of even quartics whose monomials are squares, 

£ ■= {<?40o£ 4 + q 4 y A + <?004^ 4 + 2g 2 2o£ V + 2q 202 x 2 z 2 + 2q 022 y 2 z 2 : [g 400 , • • • , ^022] e P 5 } • 

(Note the coefficients of 2 on the last three terms. This scaling tempers the coefficients 
in the equation f 5 in Figured] for the even Luroth quartics.) We show that this Newton 
polytope is a bipyramid that is affinely isomorphic to 



(»> —{(§).(§)■ a)- (§)■(!)} 

We will furthermore use the numerical interpolation method of [2] to compute the equation 
for the hypersurface in £ of even Luroth quartics. 

If i\, . . . , £5 are general linear forms on P 2 , then the quartic with equation 
(5-2) l^Wn -ii + i + i + i + i) = 

contains the ten points of pairwise intersection of the five lines defined by £i,...,£ 5 . 
Counting constants suggests that there is a 14- dimensional family of such quartics, but 
Luroth showed [20] that the set of such quartics forms a hypersurface in P 14 . 

The formula (I5.2p exhibits the Luroth hypersurface CH as the closure of the intersection 
of a general affine hyperplane Wl C C 15 with the image of the map 

g : (c 3 ) 5 — ► C 15 

(5 - 3) (4,. .,4) h-> fl^-ti 

i=l i=l 

The codimension of CH and the dimension of the general fiber (both 1) are easily 
verified using this parameterization [151 Lemma 3]. In particular, we used the method 
of pj)] described in § [2] with Bertini pQ to compute a witness set for CH = ^(C 15 ) R 971. 
This witness set verifies that the degree of CH is 54. As shown in [TB], this witness set 
also provides the ability to test membership in CH by tracking at most 54 paths. 
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The space C 15 of quartic polynomials has coordinates given by the coefficients of the 
monomials in a quartic, 

^2 QijkX l y : 'z k = q mQ x 4 + q 3W x 3 y + q 301 x 3 z + g 2 2o^V + Q2nX 2 yz 

i+j+k=A +q202X 2 z 2 + qi 30 xy 3 + qn\xy 2 z + q U2 xyz 2 + gios^z 3 

+go4o2/ 2 + qo3iy 3 z + q 02 2y 2 z 2 + g i3^ 3 + qoo^z 4 . 

In theory, we may use the algorithm of Remark [10] to determine the Newton polytope 
of CH. While difficult in practice, we may compute some vertices. For example, 

OfoitfMa (6,0,6,0,0,0,0,30,0,0,0,0,0,12,0), 

is the extreme monomial in the direction 

(3, -5, 3, 2, 3, -2, -1,4, -3, -2, 3, 1, -5, 3, -5) . 

By symmetry, this gives five other vertices, 

9400?310?1129031 , 9o40^031?211<?103 i <?040 <?130<?1 12<?301 ) 9o049oi3921l9l30 i 9o049l03?12l9310 • 

It is dramatically more feasible to compute the Newton polytope of the hypersurface of 
Liiroth quartics in the space E of even quartics. This is the face of the Liiroth polytope 
that is extreme in the direction of v, where 

v • (<?4oo, 93io, • • • , 9004) = - ^{lijk I one of i, j, k is odd} . 

Obtaining a witness set for the even Liiroth quartics, EH :— E fl CH, is straightforward; 
we reparameterize using the 2's in the definition of E and include the linear equations 

q%jk = where one of i, j, k is odd 

among the affine linear equations C: C 15 — > C 13 used for the witness set computation. 
When performing this specialization, some of the 54 points from CH coalesce. More 
precisely, six points of E% arise as the coalescence of four points each, nine points of EH 
arise the coalescence of two points each, and the remaining twelve points remain distinct. 
This implies that £% is reducible with non-reduced components. 

Numerical irreducible decomposition shows that E% consists of eight components, only 
one of which is reduced. However, as we are using witness sets for images of maps pj)] 
(as described in § [2]) the numerical computations are not performed in E, but rather on 
the smooth incidence variety of the map g (15. 3p . 

We first determine the Newton polytope of each component and then use interpola- 
tion [2] to recover the defining equation for each component. For j\, . . . , f 5 as given in 
Figure [U E7-L is defined by 

(5-4) g 4 4 00 ■ g 4 40 ■ g 4 04 • ft ■ fl ■ fl ■ fl ■ h = . 

For completeness, we used the algorithm of [16] to verify that a random element of each 
hypersurface V(/i) lies on CH. 

Observe that fi, f 2 , f 3 , and / 4 all have the same support and therefore the same Newton 
polytope, A. Every integer point of A corresponds to a monomial in these polynomials 
and all are extreme. The Newton polytope of f$ is 4A and it has 65 nonzero terms, 
which correspond to all the integer points in 4A. Thus the Newton polytope of E~H is 
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r 2 2 2 o 

Jl — g400<?040<?004 — <?400<Zo22 ~~ <?040<?202 — <?004<?220 ~ Z 9'220 , ?202 , ?022 

f2 = 9400Q'040Q , 004 ~ <7400<?(}22 + 3<?040<?202 ~~ <7004<?220 + 2q , 220Q'202Q , 022 

f% — 9400Q'040 ( ?004 + 3^400^022 — 90409202 — <?004<7220 + 2^22092029022 

fi = 940090409004 — <?400<?022 — <?040<?202 + 3<?004<?220 + 2q , 220?202Q , 022 

/ 5 = 2401g4 00 go 40 go 04 - 196g4 00 g0409004<?022 + 1029400^040^004^022 _ 4g4 00 g040?004Q , 022 

+ <?400<7022 ~ 1 96g4 00 go40?004?202 ~ 19 6gfoo9o4o9o04<?220 + 8 40g|oo9 , 04o5 , 004?220?2023022 

— 820g|ooQ , 040^0049202?022 _ 820g40 g0409004?2205 , 022 + 232g4 00 go40^0049220?202q , 022 

3 2 2 4 3 2 2 4 /13 5 

— 12g 4 o q , 04o5 , 004Q , 202%22 ~~ 12g 4 oo5 , 040%04522o5 , 022 ~ 40g 4 ooq , 040Q , 004Q , 220Q , 202Q , 022 

o 26 3 26 3 7 2424 

+ 4g 40 o<?040<?202<?022 + 4g 40 o<7oo4<?22o ( 7o22 ~~ 8g 4 o g22o<?202<?o22 + 1 02g 4 oog 4o%o4Q , 202 

— 820g4o go4o?oo4922o'?202 + 232g4\, golK)9oo4922o9202<?022 - l 2 ^^^^*?^ 

+ 102g4 00 go4o5 , 004Q , 220 + 232g4oogo4oQ , 004?220?2029022 + 128^400^040^004 5 , 220?202?022 
22 33 2244 2 342 

— 80g4o q , 040?00492209202%22 + 6g , 400%4o5 , 2029o22 — ^Q400 ( 1040 ( lo04,Q220'lo22 

2 23 3 2 224 /i 2 35 

— 80g 4 ooq , 040Q , 004922o9202'?022 + 220g 4 oog040Q , 004Q , 220?202%22 — 24g 400 <7040 ( 7220 ( 7202 ( ?022 

+ 6g400?00492209022 ~ 24g 4 \ )0 go04g220<?202<?022 + 24g 4 \ )0 g220 , ?2029o22 

— 4g 4 oog , 040 ( ?0045 , 202 — 12g400?040 ( ?004?220 ( ?202 — 40g 40 ogo40 ( ?004<Z220 ( ?202 ( ?022 

362 ^34^ 2233 

+ 4g400<?040<?202<?022 — ^^OO^O^wfeo'^ ~ 80g 4 oogo40 , ?004 , ?220 ( /202 , ?022 

+ 220g400q , 040?004g220Q , 202Q , 022 — 24g 4 oogo4oQ'220 ( ?202 ( ?022 — 4g 40 og0409o04?220 

35 2422 333 

— 40g400g040Q , 004<?2209 , 202<?022 + 220g400g040?004<?220<?202 , ?022 ~~ 272g 4 oog040<?004<?220<?202<?022 

+ 48g400Q , 0409220Q , 202?022 + 4^400^0^ <?220<?022 ~~ 24g 40 0<7o04 ^LoO^^c^ 

+ 48g400g004g2203202<?022 — 32g 4 oog22o32029o22 + 9o40<?202 + 4^040^004 <?22o9202 

— 8go 40 g220g , 2 7 02 , ?022 + 6g0409004<7220<?202 ~ 24go 40 g 04g220 , ?202 , ?022 + 24<7 040 <? 22 o3202 , 2 , 022 
+ 4g040?004'?220'?202 ~~ 24go40?0049220?202?022 + 48g 409004g220?202?022 

—32 q , 040Q , 220?202Q , 022 + 9004^220 ~~ 8<?004 Q'220?202g022 + 24^004 ?220<?202<?022 

— 32q , 004g2205 , 202%22 + 16Q , 220 < ?202<?022 ■ 

Figure 1 . Polynomials denning EH 

14A + a, where a is the exponent vector of ^oo^cLo^oW To complete the identification of 
M{S%), consider the integer points {O, A, B, C, D} of A, which are on the left in Tabled] 
Replacing {O, . . . , D} by their differences with O gives the points o, a, b, c, d on the right 
in Table [TJ Note that a + b + c = d. Projecting to the first three coordinates is an 
isomorphism of the integer span of a, b, c with Z 3 , and shows that A is affinely isomorphic 
to the bipyramid (15. ip 

Using (15. 4p . we can determine which Edge quartics [9j|23] are Liiroth quartics since the 
family of Edge quartics £T> is contained in £ with 

£V := {V(s(x 4 + y 4 + z 4 )-t(y 2 z 2 + x 2 z 2 + x 2 y 2 )) ■ [M]^ 1 }. 
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Table I. 





<?400 


<?040 


#004 


%22 


Q202 


<?220 














1 


1 


1 


A 


1 








2 








B 





1 








2 





C 








1 








2 


D 


1 


1 


1 












Vertices of A 





<?400 


<?040 


<?004 


Q022 


Q202 


Q220 























(1 


1 








1 


-1 


-1 


b 





1 





-1 


1 


-1 


c 








1 


-1 


-1 


1 


d 


1 


1 


1 


-1 


-1 


-1 



Identifying £V with P 1 , and evaluating at (15. 4p gives the equation for £V D CH, 
s 12 (s + t) 4 (2s-t) w (7s + t)(2s 2 + st + t 2 ) 6 (28s 3 + 8s 2 t + 3st 2 + t 3 ) 3 = 



Set u := V297 + 24\/159. Besides the point [0, 1], the eight points corresponding to 
Edge quartics that are Liiroth quartics are 



tl = 


-1, 


t 2 = 


2, 


h = 


-r, 


U = 








h = 






f (1+ ' 


h = 






1 ( 2 


tr = 


— (oj — 




ts = 


1 f 2 









- 15 + v 73 ^ (lu 2 + 15)) , 

- 15 - V^3 (to 2 + 15)) . 

In particular, there are four real values ti,t 2 ,t 3 ,tQ and four nonreal values t 4 ,t 5 ,t 7 ,t 8 . 
The Edge Liiroth quartic corresponding to [0, 1] has three real points, each of which is 
singular. Also, except for t = t 2 = 2, which is the union of four lines 

x — y + z = x — y — z = x + y — z = x + y + z = 0, 

the Edge Liiroth quartic corresponding to [1,U] is smooth with no real points. 



6. Conclusion 

We presented two algorithms for computing the Newton polytope of a hypersurface % 
given numerically. The first assumes that we may evaluate a polynomial defining % while 
the second uses a witness set representation of H. The second is illustrated through the 
determination of the polynomial defining the hypersurface of even Liiroth quartics (which 
gives a face of the Liiroth poytope), along with some other vertices of the Liiroth polytope. 
Implementing these algorithms remains a future project. 
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